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Abstract 

The problem of two Aharonov-Bohm (AB) vortices for the Helmholtz equation is examined in detail. It is demonstrated that 
the method proposed in [J. M. Myers, J. Math. Phys. 6, 1839 (1963)], Ref. [4], for diffraction on a slit can be generalized to get 
an explicit solution for AB vortices. Due to singular nature of AB interaction the Green function and the scattering amplitude 
for two AB vortices obey a series of partial differential equations. Coefficients entering these equations, in their turn, fulfill 
ordinary non-linear differential equations whose solutions can be obtained from a solution of the Painleve V (or III) equation. 
The asymptotics of necessary functions for very large and very small distances between two vortices are calculated explicitly. 
Taken together, it means that the problem of two AB vortices is integrable. 
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I. INTRODUCTION 


The Aharonov-Bohm (AB) effect [T] , [2] is one of the most striking distinguishes between quantum and classical 
words. In a nutshell, it states that a quantum particle feels electro-magnetic potentials even though no classical forces 
exist. Its first description can be traced to the paper of Ehrenberg and Siday [Tj, but it is only after the seminal work 
of Aharonov and Bohm [2] that this subject attracts a wide attention. The success of that paper can be attributed 
to the fact that in addition to a general discussion of the phenomenon the authors presented a clear-cut analytic 
calculation of physical scattering on one singular AB vortex thus validating common arguments. 

Today there exists a huge literature about this effect but, surprisingly, analytically results are rare. In [5] a 
diagrammatic-like series for the amplitude of scattering on a few AB vortices had been proposed but for real energy 
it is similar to a formal multiple scattering expansion and hardly can be used for calculations. 

The purpose of this paper is to investigate the problem of scattering on two AB vortices. The principal result is 
that this problem is integrable and the calculations of the Green function and the scattering amplitude can be reduced 
to a solution of a series of differential equations whose lowest level includes the Painleve V (or III) equation. The 
method used in derivation of these results is a generalization of the one proposed in Ref. [4] where the diffraction on 
a finite slit has been treated. It is based on the point-like nature of the AB potential which permits to fix solutions 
by fixing its behavior near vortex positions. It means that only a few constants uniquely determine the full solution. 
Using different transformations commuting with the Laplacian leads to a system of equations for these constants. 
Besides equations, it is necessary to know the values of different quantities at small and/or large distances between 
vortices. For large vortex separation it can be done by perturbation series and for small distances between vortices it 
is achieved by the using Riemann-Hilbert method. 

The plan of the paper is the following. Section [Tl] is devoted to a general discussion of the problem of scattering 
on AB vortices. Special attention in this Section is focused on the uniqueness of the solution and, in particular, on 
the fact that any solution obeying all boundary and radiation conditions but without in-coming incident waves is 
identically zero. In Sections Hl||V it is demonstrated how the arguments of Ref. can be generalized to the case of 
two vortices. First of all, a set of auxiliary functions independent of incident fields with prescribed singularities at 
vortex positions are introduced. These functions are analogues of the Hankel functions for one-vortex problem and 
play a prominent role in what follows. The main idea of Ref. jTj is that there exists a group of differential operators 
which commute with the Lagrangian and cancel the incident field. Transformed solution is non-zero as the action of 
these group operators change boundary conditions near vortices. But these changes can be compensated by a suitable 
linear combination of new functions. In this manner one gets a set of equations for unknown functions. Calculation 
of group commutators done in Section |IV| p ermits to find equations for all necessary functions. In Section |III| the 
Green function is discussed and in Section |V| this procedure is done for the scattering amplitude. To use the obtained 
equations one needs to find the asymptotics of correct solutions at small and/or at large separation between vortices. 
This is achieved in Section |VI| where explicit forms of the solution when the distance between vortices tends to zero 
and to infinity are obtained. Section |YlI| is a summary of the obtained results. The relation of these results to the 
theory of holonomic quantum held [2j is in short discussed here. Appendix [A] is devoted to the proof of the uniqueness 
of the solution and to derivation of the reciprocity relation for the scattering on two AB vortices. In Appendix [B] 
properties of one-vortex solution are briefly discussed. 

To diminish the paper size only the main steps of derivations are presented and details are often omitted. 


II. GENERAL CONSIDERATIONS 

The AB vector potential, A^, is a pure gauge potential, A M = d^cj), and can be removed by a gauge transformation. 
Nevertheless, the existence of AB vortices manifests in non-zero circulation along any closed contour encircling only 
vortex j 

j) A^dx^ = -27ray, (1) 

where aq is the magnetic flux associated with the vortex (we assume that a ^ integer). 

The existence of a non-zero circulation implies that after potential removing from each vortex emanates a line of 
phase discontinuity (the cut) denoted Cj such that the function and its normal derivative on the both sides of the cut 
differ by a phase 

4- + (x,0) = e^^-MuO), 0 v $+(a:,O) = e 2 ™’d y V-(x,0), x e C 0 . (2) 

Here x and y are coordinate respectively along and perpendicular to the cut. Each cut has two different sides which 
can be connected by a contour encircling one or more vortices. The cuts can be chosen arbitrarily and wave functions 
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FIG. 1. Cuts for the two-vortex problem (red lines). Black solid lines are contours encircling the vortices, 
indicates large radius contour around all vortices. 


with different cuts are gauge equivalent. In Fig. |T]a convenient choice of cuts for two AB vortices used throughout the 
paper is sketched. In this case the two cuts coincide along the negative x-axis and boundary values of wave function 
and its y-derivative on the cut are related as follows 

$+(*,0) = e 27rix(:c) ’F_(x,0), 0 v tf + (x, 0) = e 2nix ^d v ^_{x, 0), (3) 

where piece-wise constant function %(x) is 


+ « 2 , x < 0 
oq, 0 < x < L 
0, x > 0 


( 4 ) 


and subscripts (±) corresponds to limit y —> 0 from positive and negative values of y respectively. The problem of 
one vortex has been solved in [2] (and is shortly reviewed in Appendix [B|. 


As any problem of diffraction, the scattering on the AB vortices corresponds to finding a wave function ^(x) with 
the following properties: 

(a) \I/(x) is the sum of an incident wave T lnc (a;) which includes all in-coming waves and a reflected out-going wave 
4- ref (x), 

1>(x) = ^ inc (f) + T ref (f) . (5) 

The choice of the incident wave is dictated by the problem considered. When one is interested in the Green 
function, the incident wave is the Green function of the free Helmholtz equation 

V™{x) = ^H£\k\x-x'\) (6) 

with x' being the source point. For the problem of the plane wave scattering 

^ inc (f) =e ifercos(e -^, (7) 

where r and 9 are polar coordinates of point x and (j) is the angle of incidence. 

(b) The reflected field obeys the Helmholtz equation 

(d 2 x +d 2 y +k 2 W ef (x)= 0 (8) 

everywhere in the plane (x, y) except the cuts. 

(c) At the both sides of the cuts the full wave function and its normal derivative are related as in Eq. |2| (or for the 
cuts in Fig. [l] as in Eq. <§)• 

(d) At large |x | the reflected field obeys the out-going radiation condition which is legitimate to choose in the form 

lim / |c> r 'I' ref — ifcT ref | 2 ds = 0 (9) 

fi-*- 00 Jc R 

where the integration is performed over a big circle Cr of radius R which includes all vortices and s is the length 
along this circle (see dashed circle in Fig. [TJ 
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(e) Vortices are considered to be impenetrable. It means that the full wave function tends to zero at vortex positions. 
More precisely, in a small vicinity of each vortex the full wave function should have the following behavior 

\I/(x) —> a,j(x — Lj + i y) aj + bj{x — Lj — i y) 1 ~ aj ■ (10) 

x—>Lj 

As the integer part of flux does not change boundary conditions we may and will consider fluxes in the interval 
0 < aj < 1. 

From physical considerations it is clear that the solution obeying all conditions (a)-(e) is unique or, which is the same, 
if a function fulfills all these conditions with zero in-going incident field, it is identical zero. The proof of this fact 
can be done by a generalization of usual arguments developed for diffraction problems (see e.g. [S] and references 
therein). For completeness, we present in Appendix [A] a brief demonstration of uniqueness stressing the necessity of 
all requirements (a)-(e). 

The usual way of solving the problem with two vortices is to represent the reflected outgoing field as a sum of single 
and double layers along the cuts 

\F ref (cc, j/) = f H™ (ky/(x - t) 2 +y 2 ) y.{t)dt + d y f (ky/(x - t ) 2 + y 2 ) v(t)dt , (11) 

J cuts ' ' J cuts ' ' 

where h'Y (x) are the Hankel function of the first kind and zero order. 

Functions /i(f) and u(t) are piece-wise functions on the cuts which have to be determined from the boundary 
conditions ©• The calculations are standard (cf. e.g. [7]) and lead to the following system of equations 

v{x) = g tan7rx(x) / HY\k\x — t\)y(t)dt + F(x, 0) , (12) 

J cuts 

/x(x) = — | tan7rx(x) (px + k 2 ) j H^\k\x — t\)v(t)&t + d y T{x,§) , (13) 

J cuts 

where (provided that the incident wave T lnc (x) has no phase jumps) 

.F(x, 0) = \ tan7r\(x) T inc (x, 0), d y F(x, 0) = ^ tan7rx(x) 9 s $ mc (x, 0) . (14) 

This approach is well suited for numerical calculations. To progress in analytic treatment, we generalize in the next 
Section the method of Ref. U which has been developed for the diffraction on a finite slit. 


III. GREEN FUNCTION FOR TWO-VORTEX PROBLEM 


To calculate the Green function for the two-vortex problem it is necessary to fix the incident field as in Eq. (Hi- 
Due to dependence of the source coordinates, x', the asymptotics of the Green function denoted by G(x,x') is 


G(x, x') 


x—>Lj 


j(x') (x - Lj + i y) aj + bj(x') (x - Lj - i y) 


l—a-j 


(15) 


The choice of fractional power branch is dictated by the choice of the cuts. 

Let us introduce auxiliary functions Aj (x) and Bj (x) independent of the incident field which obey all the conditions 
([8])-(10l except that at one vortex indicated by j they have the following asymptotic behavior 


Aj(x) 


x—±L- 


d Lj (x - Lj + iy) a ^ = - 


B j(x) d L (x- Lj-iyY 

x^Lj 


(x — Lj + \y ) 1 a i 

1 - OLi 


(16) 

(17) 


(x — Lj + i y) a i 

These functions are uniquely fixed by these conditions and will play an important role below. 

The uniqueness of the solution together with point-like character of boundary conditions (15) permit to find different 
equations between functions aj , bj , Aj , Bj j3] . Indeed, assume that there exists an infinitesimal transformation S 
commuting with the Laplace operator such that the incident field is invariant, d4' mc (x) = 0. Then the function 
i5T(x) is a solution with zero incident field. Nevertheless, such transformed solution is not zero because, in general, 
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the symmetry transformation 6 changes the behavior near one or many vortices. But this change can be compensated 
by a suitable chosen linear combinations of functions Aj and Bj . It means that the function 


\- c i A o (%) + d .i B .i )] 


(18) 


obeys all conditions (|8j)-( 10) but with zero internal field. By uniqueness it has to be identically zero, 

<5\l/(:r) + y^[cj-Aj(aT) + djBj(x)\ = 0 . 


(19) 


Combining together different transformations leads to sufficient number of equations which permits in the end to 
reconstruct the full solution. 

When we are interested in the Green function for a system of vortices, the incident field ([6]) has the following main 
symmetries 

• change of vortex positions: S = <9r., 

• translational invariance: 5 = d x + d x >, 5 = d y + d y i. 

• rotational invariance: 6 = = xd y — yd x — Ljd y . 

In the subsequent Sections these transformations and their combinations are considered and equations for the Green 
function of two AB vortices are derived. 


A. Derivatives over vortex positions 


From (15) and definitions (161 and it follows that the combination 

dLjG(x,x') — a,j(x')Aj(x) — bj(x')Bj(x) 

is zero at all vortices and, as was discussed above, is identically zero. 

Therefore derivatives of the Green function over vortex positions are 

dLjG(x,x') = cij(x')Aj(x ) + bj(x')Bj(x) . 

Assume for a moment that 0 < a < 1/2 then dj(x') and bj(x') can be calculated from the limits 

aj(x') = lim (x — Lj)~ aj G(x,x'), bj(x') = lim (x — Lj) aj ~ l [G(x,x') — a,j{x'){x — Ai) Qj ] 

x—^Lj x—>Lj 

Differentiating this limit over L 2 for j = 1 and over L\ for j = 2 one concludes that 


9l 2 


ai(x') 

bi(x') 


Pi ei 

Si Ci 


a 2 (x') 

b 2 {x') 


<9z,i 


a 2 (x') 

b 2 {x') 


Pi e 2 
S 2 C 2 


ai(x') 

bi C x ') 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


where all matrix elements are independent on space coordinates. 

Calculating the mixed derivatives d1 lL2 G(x,x') one gets that the derivatives of A,; and Bi are expressed through 
the same constants as follows 


9l 2 


Ai(x) 
Bi (x) 


p 2 5 2 

e 2 C2 


A 2 (x) 

B 2 (x) 


9l 2 


A 2 (x) 

B 2 (x) 


Pi Si 

ei Ci 


Ai(x) \ 
Bi{x) J 


(24) 


B. Translational invariance 

There exists a few simple consequences of translational invariance. First, constants like pj, €j , 5j, Q depend only on 
the difference L = L\ — L 2 . Second, functions depended only on x or x 1 like fi = cii,bi, Ai,Bi obey the equation 

{d Ll + d L2 + d x ) fi(x) = 0 . (25) 
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This equation permits to calculate derivatives dLifi in contrast to Eqs. (23) which determine derivatives 8l ft only 
with j ^ i. 

Consider now the change of two coordinates simultaneously. Comparing the behavior near two vortices (15) one 
concludes that 


\(d y + d v >) G(x , x') = ai(x')Ai(x) - b 1 (x')B 1 (x) + a 2 (x')A 2 (x) - b 2 (x')B 2 (x) . 


Differentiating this expression by L x and using Eqs. (24) and (25) leads to the identity valid for all x and x 1 


A 1 {x){d x i + \d y ‘)ai(x') + a 1 {x'){d x + id v )A 1 (x) - B 1 (x)(d x > - idy^b^x') - h(x')(d x - i d y )B 1 {x) = 
- Ai(x)[fi x a 2 (x') + eife 2 (x')] - ai (x')\/3 2 A 2 (x) + 6 2 B 2 (x)\ + A 2 (x )[/3 2 ai(£') + e 2 bi(x')} 

+ a 2 (x')[j3\Ai(x) + Sib^x')} - B 2 (x)[5 2 a x (x') + ( 2 bi(x')] - b 2 (x ')[eiAi(£) + CiB x (x )] . 


Substituting in this expression the most general linear relations between derivatives and unknown functions 


( d x + id y )A 1 (x) = giA x (x) + g 2 B 1 (x) + h 1 A 2 (x) + h 2 B 2 (x) , 

(■ d x - id y )B x (x) = g 3 A 1 (x) + g i B 1 (x) + h 3 A 2 (x) + h 4 B 2 (x) , 

(d x i + idyi) a x (x') = G\ai{x') + G 2 b 3 (x') + Hia x (x') + H 2 b 2 (x') , 
(d x > - idyi) h(x') = G 3 ai(x') + G 4 b x (x') + H 3 a x (x ') + H 4 b 2 (x') , 
and collecting identical terms proves the following formulae 


(d x + id y )A 1 (x) = g 1 A 1 (x) + g 2 Bi(x) - 2 5 2 B 2 (x) , 

( d x - id y )B x (x ) = g 3 Ai (x) + g 4 B x (x) - 2e 2 A- 2 (x) , 

(d x i + idyi) ax(x') = -gi a 4 (x') + g 3 b x (x') - 2e 1 b 2 (x') , 

(d x i - i d v ')bi(x') = g 2 ai(x') - g 4 h(x') - 2<5ici 2 (£ / ) . 

Similarly, differentiating Eq. ( |26| ) by L 2 gives rise to the relations 

(d x + id y ) A 2 (x) = fiA 2 (x) + f 2 B 2 (x ) - 2S 1 B 1 (x) , 

(d x - idy) B 2 (x) = f 3 A 2 (x) + f 4 B 2 (x) - 2<nA x (x) , 

(■ d x > + id y i) a 2 (x') = -/i a 2 (x') + f 3 b 2 (x') - 2e 2 b 1 (x’) , 

(d x > - id y i)b 2 (x') = f 2 a 2 (£') + f 4 b 2 (x’) - 2S 2 a 4 (£') • 

Here gj and fj with j = 1, 2, 3,4 are undetermined constants depended only on L = L x — L 2 . 


(26) 


(27) 


(28) 


(29) 


(30) 


C. Rotational invariance 


The mutual rotation of x and x' around any vortex is evidently the invariance transformation for the incident field 
Q. Such rotations are generated by operators 



L ( g j) = xd y - yd x - Ljdy , L^) = x'd v > - y'd x > - Ljd y i . 

(31) 

One has 




(4 J) + l$>) H™ (k^(x-x'Y + (y-y') 2 ) = o • 

(32) 

Exactly as 

it has been done above one concludes that 



i G(x,x ') = Lai(x')A x (x) — Lb x (x')B x (x) , 

i + LgG(x , x') = L a 2 (x')A 2 (x) — L b 2 (x')B 2 (x) 

(33) 


with L = L i — L 2 . 

Differentiating these equations by Lj one finds the following relations (similar equations for functions ai and bj are 
not presented) 





, ( 34 ) 


where Mj are 2x2 matrices with undetermined coefficients. 
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IV. MEANING AND DETERMINATION OF COEFFICIENTS 


Eqs. (29), (30) and (34) are concise consequences of symmetries of the incident field <§• But these equations 
contain many coefficients depended, in general, on the distance between vortices, L. 

The first simplification comes from the fact that functions aj and bj for two vortices can be expressed through 
functions Aj and Bj by the reciprocity relation (A14) proven in Appendix [A] It states 


G(x, x') = G{Sx ', Sx), S(x, y) = ( x , —y) . 


Using Eqs. (21) and (33) and taking into account that F(Sx) = — SLg J> F(x ) one concludes that 


■ U) 


Ai(x) 

Bi{x) 


0 ii 

h 0 


ai(Sx) 

bi(Sx) 


A 2 {x) 

B 2 {x) 


0 t 2 
t 2 0 


a 2 (Sx) 
b 2 (Sx ) 


Substituting these values to Eqs. (24) and comparing with Eqs. (23) gives 

= P 


( P 2 £2 

V ^2 C 2 


Cl £ 1 

Si Pi 


p = 


(35) 


(36) 


(37) 


These relations mean that one can consider equations only for functions Aj and Bj. Values of constants tj and p are 
calculated in Appendix (|B|). 

In the Myers method Hj used in the previous Sections constants appeared in different formulae are just the most 


general coefficients of expansion (cf. e.g. Eqs. (28 1 ). They also can be interpreted as sub-leading terms of expansion 
of auxiliary functions Aj (x ) and Bj (x ) in smaf 


/ M{x) \ 

( 

A 2 (x) 


Bx(x) 

_ 

x — yL\ 

\B 2 (x)J 



1 vicinity of vortex positions 


Pi(x - Lx + i y) ai +5i(x- Li- iy) 1_ai 
d Ll (x — Lx — iy) 1- " 1 - y 3 (x - L ± + i y) ai 
\ ex(x - Lx + iy) ai + (x(x - Lx - iy) 1 - 011 


(38) 


and 


I A x{x) \ 

( 

A 2 (x) 


Bx{x) 

_ 

X—^L<2 

\ B 2 (x) ) 

\ 


d L2 {x - L 2 + i y) a2 - ^f 2 (x - L 2 - iy ) 1_ “ 2 
ex(x - L 2 + i y) a2 + &{x - L 2 - iy ) 1-a2 
\ 8l 2 {x - L 2 - iy) 1_Ql - 5/3(2; - L 2 + iy )“ 2 ) 


(39) 


Postulating these relations and comparing the dominant singularities at the vortices one can prove all Eqs. (23), (24), 
(29), (30), (34) as it has been done in 0. 


The first series of relations between these coefficients is obtained by differentiating the both sides of Eqs. (29) and 
(30) by Lj. As all derivatives are known by previous formulae one gets many interrelations between the constants. In 
particular, it follows that 


8l9\ = d L g4 = d L fi = d L / 4 = 0 . 


(40) 


It means that constants 51 , 54 , fi,f± are independent on L. From the asymptotic behavior of two vortices at large 
separation given in Appendix [B] one gets that 51 = 54 = /1 = / 4 = 0. We rewrite the above equations with these 
values 


(■ d x +id y ) 


a x(x') 
a 2 (x') 


and 


(d x +id y ) 


Ax(x) 

A 2 (x) 


53 -2d \ / bx{x') \ 

- 2 e 2 h ) \ ^ 2 {x') ) 


{d x -i d y ) 


bx{x') 

b 2 (x') 


52 —25i 

— 2^2 f 2 


ax(x') \ 
a 2 \x’) ) 


, (41) 


52 -262 \ ( Bx{x) \ 

-2 Sx h ) V B a(^) ) ’ 


(d x - i d v ) 


Bx{x) 

B 2 (x) 


53 —2e 2 

-2ei h 


Ax(x) \ 

A 2 (x)J 


• (42) 


Applying (d x ± id y ) operator to these equations and using the fact that each function f(x) has to obey the Helmholtz 
equation, {d x + d y + k 2 )f(x) = 0 , one finds that 

5253 + 4e 2 <5i + fc 2 = 0, £ 2/2 + ^253 = 0 , S 2 f 3 + £2 52 = 0, £ 1^2 = < 5 i£ 2 - (43) 
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Differentiating (421 by Lj gives 


& = 2(Cife + ft*), 53 = 2(e 2 /3i + £ 1 (^ 2 ), /2 = -2(<5i £ 2 + 82 P 1 ), fe = — 2 (£ie 2 + /3 2 ei) 

2^2 = / 2 /?2 - 52 C 2 , 2 e 2 = / 3 £ 2 — 5 3 / 3 2 , 2^i =/ 2 Ci — 52 / 3 i, 2 ei = f 3 pi — 5 3 Ci- (44) 

Here and below the dot indicates the derivative over L. 

In a similar way the differentiation by Lj the rotation equations (34) gets many other relations. In particular, the 
both matrices Mi and M 2 are constant. From the solution of well-separated vortices (see Appendix [B]) it follows that 
the both matrices are diagonal. For simplicity we impose this condition from now. Therefore the rotations have the 
form 


i Lg ^A\ — mi-Ai — L/3 2 A 2 + LS 2 B 2 , i Lg ^B\ — niBi — Le 2 A 2 + LfeB 2 

i Lg ^A 2 = m 2 A 2 + L(3\Ai — L8\B\ , iLg ^B 2 = n 2 B 2 + Le±Ai — LfeBi 

In the indicated way one gets the following list of equations 

(Lei)' = {mi - n 2 )ei - Lfeg 3 , (Le 2 )' = (m 2 - ni)e 2 + Lfefe, 

(L(i)' = (n 2 - m)Ci - Lei 5 2 , (AC 2 )' = (ni - n 2 )fe + Le 2 / 2 , 

(L/3i)' = (mi - m 2 )/?i - L6ig 3 , (L/3 2 )‘ = (m 2 - mi)/3 2 + LS 2 f 3 , 
(LJi)' = (m 2 - ni)<5i - L/3ig 2 , (L<5 2 )' = (mi - n 2 )S 2 + L/3 2 fe. 

According to Eqs. (B12) of Appendix [B] the values of mi and are 

mi = 1 — ai, m = — ai, m 2 = l — a 2 , n 2 = —a 2 . 


(45) 


(46) 


(47) 


As was discussed above one has to take into account only equations related with one vortex (say Li). We rewrite 
them again 


<7 — wei — Ci53> 


Cl =- £i52, Pi = ^^i - *103, 


di — — j-Si — Pig 2 (48) 


with 7 = a 2 — ai. 

To these equations one should add the following ones 


27 , 


5253 + 4e 2 <5i + k — 0, Pie 2 g 2 ~ p 2 ^ig 3 + ^S 2 e 1 — 0. 


(49) 


The second of these equations is a consequence of Eqs. (46) when relations (37) are imposed. 
It is convenient to introduce the following notations 


V = S 2 ei = Sie 2 , v = f 3 6 2 = — 5 2 e 2 , 
Direct check proves that these variables obey the equations 


w = g 3 Si = -fee 1 . 


(50) 


y = piv~p 2 w, Pi = 1 Pi -w, p 2 = - 7 + 1 p 2 +v, w = -jW + Pi(k 2 + 8 y), v = j-v - p 2 (k 2 + 85 ) (51) 

with 7 = mi — m 2 = a 2 — ai . 

Eqs. (49) are equivalent to the existence of two integrals of motion 


27 

vw — y(^y + k 2 ) = 0 , Piv + p 2 w - —y = 0 . 


(52) 


Introducing new variable, 2 = P 1 P 2 , one gets the following equations 




7 1 

Piv = -y+ -y, 


o 7 1 • 

P 2 W = -y- -y. 


(53) 


Combining with Eqs. (52) leads to the final system of equations 

2 . 1 


2 = --2 + 5 , 


4 V - j^y + zy(Ay + k ) = 0 . 


(54) 













Substitutions 


k 2 

y = - t f ’ 


Z- 7 2 

4L 2 ’ 


L = I 


transform Eqs. (541 to the form 


(55) 


—Y' 2 - 7 2 F = ZY(Y - 1), 


Z’ = x 2 Y’. 


Removing Z gives one non-linear equation for Y 


Y" = [ -L 


1 


Y' 


2 j 2 Y 


k 2Y 2 (Y- 1) 
Finally one more change of variables 

Y = 


Y' 2 - — + 2Y(Y-1)- 

x x z (Y — 1) 


V 


= Vt 


V-V 

transforms this equation to the canonical form of Painleve V equation (see e.g. [14] ) 


d 2 E 

f 1 i M 

(dVy 

2 IdE 

d t 2 " : 

^2E ' V-lJ 

V d t) 

t d t 


2 1 2 


General Painleve V equation m is 


y" = y ' 2 


l 


l 


2 y y- 1 


y' (y- !) 2 r 

t f 2 


ay - 


y . Sy(y + 1) 
't + ^T 


Therefore Eq. (|59|) is the Painleve V equation with the parameters 

0 = 0 , 


« = b 2 ’ 


e — 2 , 


5 = 0. 


(56) 


(57) 


(58) 


(59) 


(60) 


(61) 


The Painleve V equation with <5 = 0 can be reduced to the Painleve III equation m- Consider the following 
Backlund-type transformation for two functions Y and W 


■ _ 2 7 2Y(Y - !) 
x W '' 


W = 1 - 2Y + 1 + 27 IE + IE 2 


(62) 


Finding IE from the first equation and substituting it into the second one leads to Eq. (571 for Y. Calculating Y from 
the second equation and putting it to the first one gives 


IE 2 IE , 

IE = --+ IE 3 -T ■ 

IE x 

Standard form of the Painleve III equation is (see e.g. na 

•2 


2(l+ 7 )IE 2 -2 7 1 

x W 


r y . 3 0y 2 + e 6 

y = -h ay H-h - 

y x x y 


Therefore Eq. (63) is the Painleve III equation with parameters 

0 = 2(l + 7 ), 


a = 1, 


e = -2 7 , 


S = - 1. 


(63) 


(64) 


(65) 


The knowledge of the solution y = y(L) permits to find all other quantities by simple integration. First, z(L) is 
determined directly from the second of Eqs. (541. From Eqs. (511 it follows that 


u = \Jy(4y + k 2 ) e s , w = \Jy{Yy + k 2 ) e s , 
„ l _ -IV/L + V/2 - s a _ -tv/L-m ,-s 


\J y(4y + k 2 ) 


y/ y(4y + k 2 ) 


( 66 ) 

(67) 
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where function S = S(L) is calculated from the equation 


S = 


4q y 


L(k 2 + 4 y) 


Using Eqs. (48) one demonstrates that 


ei = C^/y, <5i = ^Jy, e 2 = pC^/y, S 2 = ^y 


where C = C(a 1 , a 2 ) is a constant which can be calculated from limiting values found in Section VI 

“ 1+ “ 2_1 r(2 - q 2 ) / sin7ra2 c - ffi (3a 1+ a a v 2 
T(1 + a.\) V sin7rai 


C(a 1 ,a 2 ) = ( - 


Other quantities can be calculated from the definition and reciprocity relation (371 

92 = 


1 


Cl = -02, C2 = P01, 

p 


V 

£2 ’ 


9 3 = 


w 

V 


( 68 ) 


(69) 


(70) 


(71) 


V. SCATTERING AMPLITUDE 


A typical scattering problem consists in the determination of wave function, 'I' scat (x). when the incident field is 
chosen as plane wave (J7|. The main physical quantity of interest is the scattering amplitude, F(9, (f>), obtained from 
the asymptotic behavior of this function 


\E- scat (f) —-> e 

\x \—>00 


i kr cos { 6 —<p) 


7 ri kr 


e Ifer ?(0, 0). 


(72) 


Here x = (r cos 9, r sin 9) and angle (f> determines the direction of the incident plane wave. 

Let us consider the limit of the Green function, G{x,x'), with x' = (Rcoscp. Rsm®) and R —> 00 . As 


H^\k\x-x'\) 


lim 


_ —ikr cos 

1 


(73) 


the scattering wave function with asymptotics (721 can be extracted from the Green function as follows (cf. [1]) 

G(x, —x') 


vp s 


(x) = lim - 

R—t 00 1 


hH { v\k\x'\ 


(74) 


Functions Aj{x) and Bj{x) defined in the previous Section (cf. Eqs. (16) and @) obey the Helmholtz equation and 
due to the radiation condition § have the following asymptotic behavior 


Aj(x) 


| a; | —Kx) V 7ri kr 


r FA9)e~ 


B j (x) 


|rc|->oo V nikr 


r G,(d)e- 


(75) 


with certain functions Fj(6) and Gj(6). 

Because of translational invariance (<9^ +9 l 2 +d x )Aj = 0, (d^ +9 l 2 +d x )Bj = 0, functions Fj(6) = Fj(6,L) and 
Gj(9) = Gj(9, L) depend only on the distance between vortices, L = L\ — L 2 (and angle 9). They have the meaning 
of asymptotics ( |75[ ) when center of polar coordinates is chosen at vortex j. To simplify notations, the arguments of 
Fj and Gj are dropped when it will not lead to confusion. 

For the scattering on AB vortices the exact wave function ^(ir) tends to zero at the vortex positions (cf. |To|)) 


T scat (£) 


x—>L■ 


CLj(x — Lj + i y) aj + bj{x — Lj — i y) 


1 —a. 


(76) 


It is plain that the operator 5 = d x — ik cos (j> gives zero when acting on the incident plane wave in Eq. (72). Therefore 

(77) 


(54' scat (x) corresponds to the zero incident field. Comparing the behavior near the vortices one gets that 

( d x — ik cos 0)\l/ scat (x ) + ^ | djAj (, x) + bjBj (, x ) 

i 


= 0 . 
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Applying this relation to asymptotic expression |72| ) one concludes that 

ifc(cos# — cos </>) = — ^ \a,jFj(0) A bjGj(O) 


i k cos QLj 


(78) 


Such type of expressions is called embedding formulae in the theory of diffraction 01 - 

Values of a,j and bj can be calculated from the limit (74) together with {77]) and the reciprocity conditions (]36|) 


a,j = lim 


dji-x') 


4i 


R ^°° ^H^\k\x'\) tj 


= -Gj (t r - (j))e ik cos + L i, bj = lim 


bj(-x') 


4i 


R ^°° ^ 


= -Fj(TT-(l))e ikcos ^ L G (79) 


For clarity the argument of Gj and Fj functions here is written as 7r — (j). In general the choice of the branch has to 
be consistent with the position of the cut. 

Finally one obtains that the AB scattering amplitude is 


nw) = - 


l r 


fc(cos 9 — cos ft) tj . 


E Gj (n - (f>)Fj{6) A Fj(n - ft)Gj{9) 


kLj (cos <p —cos 0) 


(80) 


Gj are linear combinations of Fj 


The relations derived in previous Sections induce equations for Fj and Gj. From Eqs. (42) and it follows that functions 

/ p, \ / 

= Vi 


Gr 

G 2 


Ft 

F 2 


F t 

F 2 


= V 2 


Gi 

G 2 


where 


Fl ifc l -2eie- ifcicose 


—2e 2 e ikLcose 

fs 


V 2 = - 


-i9 


ifc 


9 2 

_ 2 6 0 —ikL cos 6 


—2<5 2 e ifci cos6> 

/2 


Conditions (43) imply that ViV 2 = 1. 


Eqs. (24) and (451 signify that derivatives of Fj and Gj over L and 6 obey the equations 

/ p. \ ... „ / \ , 

d L 


F t 

G-| 


= —e 


i/cL cos 0 


and 


i«9e 

i«9e» 


Fi 

Gi 

A 2 

g 2 


C fa s 2 

V e 2 C2 


mi 0 

0 n\ 

m 2 0 
0 n 2 


F 2 

g 2 


Fi 

G 1 
F 2 

g 2 


dr 

A Fe 

A Fe 


F 2 

g 2 


= e 


—i kL cos 6 


fti tfi 

e i Ci 


Fi 

Gr 


ikL cos 9 I @2 &2 

—e 2 C2 

— ikL cos 9 ( ft 1 ^1 

£1 -Cl 


A 2 

g 2 


A-i 

Gi 


Expressing Gj through Fj by Eqs. (81) one finds 


F-| 

F 2 


= M 


Fr 

F 2 


id e 


Fr 

F 2 


= N 


F\ 

F 2 


where matrices M and N are 


M = 


IV = 


/ r>- elS 




p 2 -i„ T 


ft - 


-. — ikL cos 6 


AkL cos 6 


e' e 

2ly T J 


mi A 2iy— 


—F 


V 


“I 


@2 + 


~i/cL cos 0 


/3i a he— 


-.—ikL cos 0 


m 2 - 2iy— j 


Here the notations are the same as in Eqs. (50). 
The compatibility condition 


d L N - i d e M = MN - NM 


(81) 


(82) 


(83) 

(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 


(89) 


is equivalent to Eqs. (51) which is another way of their derivation. 
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VI. SOLUTION AT SMALL AND LARGE VORTEX SEPARATION 


To really use equations derived in the previous Sections it is necessary to know the values of all variables at a certain 
point. In this Section it is demonstrated how to find wave function and scattering amplitude for two AB vortices 
when the separation between them, L , is small and large with respect to the wavelength. The cases of small-vortex 
separation with opposite fluxes (i.e. o ?2 = 1 — ai) and with two arbitrary fluxes ai and 012 require different arguments 
and are discussed separately. 


A. Two vortices with opposite fluxes at small distances 


The vortices with opposite fluxes is considered first. Let the vortex with flux a be in the point L\ = L and the 
second vortex with opposite flux be at L 2 = 0 (0 < a < 1). 

The full wave function is represented as in Eqs. (11) 

^{x,y) = ^inc(x,y) + J H™ (*V (x - t) 2 + y 2 ) y(t)dt + d v J (ky/(x - t) 2 + y 2 ) v{t)dt . (90) 

Equations (121 and (13) in this case take the form 

'Linc(z, 0) + / H^\k\x - t|)/z(t)dt 
Jo 


v(x) = 

M(*) = 


tan 7 ra 

2 

tan 7 ra 


inc(^5 0) 




H^\k\x — t\)v(t)dt 


(91) 

(92) 


\da ; 2 

The main simplification for small distance vortices comes from the fact that when condition 

kL < 1 (93) 

is fulfilled one can substitute in the above equations the asymptotics of the Hankel function at small arguments [B] 

2 i x 


Hq 1 ' 1 (z) -a — fin ^ + y) + 1 + 0(x 2 In a;) 
7 r V 2 / 


(94) 


where 7 = — 4/(1) is the Euler constant and drop the term proportional to k 2 in Eq. (92). 
After these approximations the equations for v{x) and y(x) become 


v(x) = 
y(x) = 


tan 7 ra 

2 

tan 7 rci! 


f^inc (x, 0) T 


2 i 


In 


dy'&incix, 0 ) - 


2 i d 2 
7 T da : 2 


| 

L 

L 

In 


y{t)dt + 
x — 1\ 


, kL 

ln -y + 7 ) + 1 


/j,(t)dt 


v(t)dt 


Deriving Eq. (95) on x, integrating by part Eq. (96), and taking into account that 

i/(0) = 0, v(L) = 0 

one transforms the above equations into the following system of equations 


v'(x) = 


tan 7 ra 


inc{x 1 0 ) 


tan 7 ra f L y,{t) 


+ 


-dt , y,(x) = 


tan 7 ra 


2 7TI J 0 t — X 

The equations are decoupled by introducing new variables 

C±(x) = v'(x) ±i n(x). 

It leads to 


^x^f f inc(^': 0 ) 


tan 7 ra f u' (t) 
x 


L —— dt 

To f - 


C+{x) = t&mra f+(x) — 


tan 7 ra f L £+(t) 


-+ 


7 r J 0 t — x 


dt , C-( x ) = tan 7 ra/_ (a;) + 


tan 7 ra f L £_(t) 


-+ 


7 r J 0 t — x 


dt , 


(95) 

(96) 

(97) 

(98) 

(99) 

( 100 ) 
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where 


f±(x) = -(d x ±id y )$f inc (x,0). 


( 101 ) 


These equations can be solved by the Rieman-Hilbert method (see e.g. [TTjb Let us introduce the following functions 
of complex argument 2 


= T=7 


c±(t) 


2tt[ / n t — Z 


d t. 


( 102 ) 


It is plain that 


1 f L C±{t), 

nij 0 t - x 


d t = $± p (x) + $^ own (x), C±(*) = $± (®) - $±(x) 


(103) 


where < t>!^ p (x) and $^ ; own (x) are the limiting values of functions (102) from, respectively, positive and negative y. 
After a little algebra one gets that Eqs. (llOOh are equivalent to 


$± p (x) = e T27ria ^f wn (x)+e^ ai smnaf ± (x) 
whose general solutions are (here \z\ > L) [IT] 


M*) = 
M*) = 


c, 


sm 7r a 


a 4 . 1 —a. 


(z — L^z 1 “ 27ri(^ — L) a 2 1 “ J Q 

C 2 sin 7TQ! 

(z — L) 1 ~ a z a 2ni(z — L) 1 ~ a z a J 0 


f+(t)(L — t) a t 

t — z 

L f-(t)(L — 
t — z 


d t , 


■dt 


with arbitrary constants C± and C 2 . 

When condition (931 is valid, functions f± (x) can be approximated by their values at 0 and 


= 

M*) = 


Cl /+(0) 

(z - L) a z 1_a 2i 

C 2 /-(0) 


(z - L) 1 ~ a z° 


2i 


1 + 
1 + 


aL — z 
(z- L^z 1 - 0 
(1 — a)L — z 
(z - L) 1 ~ a z a _ 


The same expression s can be obtained directly from (104) by imposing that $±(z) = 0(z 1 ). 
According to Eq. (103) £±(x) = 4>± p (x) — < h^ own (x), therefore when 0 < x < L 


z/(x) = - 
li(x) = - 


i sin naCi isinTraC^ sin ira /+ (0) [aL — x] sin7ra /_(0)[(1 — a)L — x] 


{L-x) a x l - a (L~x) 1 - a x a 2(L — x)“x 1_ “ 2{L - x) 1 ~ a x a 

sin 7raCi sin7r aCi i sin ira /+ (0) [aL — x] isin7ra/_(0)[(l — a)L — x\ 


(L — x)“x 1_a (L — x) 1 ~ a x° 


2 (L — x)“x 1_ 


2 {L-x) 1 ~ a x° 


Constants Cj have to be determined from conditions (97). First, it is necessary that 


f i/(x)dx = 0. 
Jo 


Second, from Eq. (95) calculated at x = 0 it follows that 

r L 


^inc(0,0 ) + ~J l n n{i)dt+ ^ J M*) di = 0. 

The Euler integral (see e.g. 1.5 of 0 ) and its derivative 

[ t x ~ 1 {l-ty- 1 dt= [ \ntt x -\l-t) y - 1 dt=[y{x)-*(x + y)] r ^ r ( y ' > 

Jo 1 (x + y) J o 


r(x + y) 1 Jo 
permit to calculate all necessary integrals analytically. 


T(x + ?/) 


(104) 


(105) 


(106) 


(107) 


(108) 


(109) 


( 110 ) 
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From (108) it follows that C '2 = — C\ and then Eq. ( |109 ) gives 

Ci = ['Fi nc (0,0) + (1 — a)f + (0)L + a/_(0)-E]A, A = 
with 7 = —'h(l) and 


1 

2w + 4i[ln(fcL/4) + 7 + /3(a) 


/3(a) = In 2 + ^4>(1 - a) + ^\H(a) + 7 . 


(Ill) 


( 112 ) 


The knowledge of v(t ) and v(t) permits to reconstruct the full wave function. In particular the scattering amplitude 
at small kL is large only for the s-wave scattering and 

l~L r n; / IT \ 1 -1 

F = / y,{t)dt = —2 ttCi = — 


H-In —+7 + /3(a) 

7T V 4 


(113) 


In il2| it has been obtained in a different manner that 

1 /1 — a\ 1 /a' 


1 /1 — a\ 1 /a\ 

/3(a) =hi4+-®(—)+-*(-)+7 + 


7 r 

2 sin 7 ra 

Using Legendre’s duplication formula (see e.g. [5] 1.2.15) it is straightforward to check that 


4- 


-) + ) + -A— = * (1 - a) + T(a) - 21 n 2 . 

V 2 / V 2 / sin 7 ra 


(114) 


(115) 


V 2 7 sin 7 ra 

Therefore these two results are identical. 

The same formulae permit to calculate limiting values of other quantities discussed in the preceding Sections. First 
one has (z = x + it/, z = x — iy) 

(d x + id y )dJ(x,y) = (d x +id y )di iac (x,y) +- f ^^-dt = [d x + i5j / )T inc (a;, y) - 4i$+(£) , (116) 

Tl J Q z-t 

(d x - id v )df(x,y) = (9a; - id v )dii nc (x,y) - - [ —^-dt = (<9 X - i9 y )4' inc (a;, y) + 4i$_(z) . 

7 T Jo Z — t 

From these expressions it follows (as it has been checked that T (£_,-) = 0) that 

^(x,y) — a 1 (x-L 1 + iy) a + b ± (x - L 1 - iy) 1_a , d>(x,y) —> a 2 (x - L 2 + iy) 1_ “ + b 2 (x - L 2 - iy)“ (117) 

x—>Li x—>L 2 


with 


a 1 = 


a 2 =e‘ 


2iC\ 




61 =- 


2 iC 1 
1 — a 


-/+( 0 )£ C 0 ” 1 , 


(118) 


2iC7i 
1 — a 


/_(0)T T“ _1 , 


b 2 = — e 1 ’ 


2iCi 


a 


+ /+(0)U 


L" 


For the Green function one has 

1 


^inc(0,0) = - L 2 1), (9, ±i9 y )4- inc (0,0) = - L 2 |)e ±i0, 


Differentiating these expressions over Lj and using definitions (|23| one finds that 


P 2 C 2 

b 2 C2 ) L-yO 


( 


a 


1 + 


2iA 


1 — a 
—2iAL _1 


■ 2a—2 


2iA L 


-i\ 


-(1 - a) 


1 + 


2iA 


a 


r —2a 


(119) 


( 120 ) 


and 


Pi ei 

^1 Ci ) L-y 0 


V 


-(!-«) 
-2iA L " 1 


2iA 


L 


-2a 


2iAL _1 \ 


1 


2iA 
1 — a 


'2a-2 


Limiting behaviors of coefficients gj and fj follow from Eqs. (|29j) and (30) 
4ia 


9 2 


L-y 0 1 — a 


-A L 


2a—2 


9 3 




L-y 0 a 


4 i( 1 - Q) AL—, A 


4ia 


L->o a 


L-y 0 1 — a 


AL 


2ck—2 


( 121 ) 


( 122 ) 


14 













































B. Two vortices with arbitrary fluxes at small distances 


General case consists of two vortices with fluxes aq and a 2 (0 < a 3 < 1) separated by a distance, L, obeying (93). 
The principal difference with opposite flux vortices (i.e. with a 2 = 1 — aq) considered above is the existence of an 
additional cut going from infinity to the vortex positions. For convenience we choose the both cuts along the cc-axis as 
in Fig. refcontour such that the function x{ x ) is as i n Eq. Q. The reflected field is chosen as in Eq. © which leads 
to Eqs. (12) and (13) for unknown functions v[x ) and /i( x ). As a consequence, one has to know these functions along 
the whole negative x-axis and not only at short cut between two vortices. To take into account the condition (93) 
explicitly it is convenient to look for the wave function of this problem in the form slightly different from Eq. (11), 
namely 

) = Vpix) + f H™ (ky/{ x - t) 2 + y 2 ) n(t)dt + d y f (ky/ (x - t ) 2 + y 2 \ v(t)dt . (123) 

Here ^(a?) is the one-vortex solution generated by the desired incident field T mc (^) (see Appendix [b]) multiplied by 
e I7rai for the cuts as in Fig. |T] It corresponds to one vortex with flux equals the total flux of two vortices 


a _ r , , _ / 0, 0 < aq + a 2 < 1 

0 - {aq + a 2 } - aq + a 2 - »7, V - j ^ x < ai + a2 < 2 


(124) 


situated at point L 2 = 0. 


Functions v(x) and n(x) have to fulfill equations (12) and (131 which we rewrite below for the convenience 

pL 

v(x) = | tan7rx(a;) / H^\k\x — t\)/x(t)dt + T(x, 0) , 

J — OO 

pL 

n(x) = — \ tan 7 rx(x) (<9 2 + k 2 ) / H^\k\x — t\)v(t)dt + d y T(x, 0 ) . 

J —OO 


(125) 

(126) 


As function ^p(x) obeys the correct boundary conditions at the cut (—oo,0], functions J"(a;,0) and d y J r (x, 0) have 
the form 


F(x, 0 ) = 0 (,) ta T 1 ^(x,O), d y J r (x, 0 ) = 0 (x) tan 71-01 d v ^fp(x, 0 ) 


(127) 


where 0 (x) = 0 for x < 0 and 0 (a;) = 1 for x > 0 . 

As the vortex separation is assumed to be small (cf. Eq.(93)), functions v(x) and n{x) in (1231 should decrease 
quickly from vortex positions so that all integrals are dominated by a vicinity of the origin. 

In such conditions one can (i) approximate the above equations usin g (|94[ ) as it has been done in the previous 
Section and (ii) use the small-a? asymptotics of function 'f'p(x) given by (|10[) 


^p( x )z—> cl(x + iy ) 13 + b(x — iy ) 1 13 


x —>0 


(128) 


with certain (known) quantities a and b (fixed by the quantity considered). For the Green functions this expansion is 
given by Eq. 


i fc^e i7rai 


i Ad-V ™ 1 


- = a ^') = - 2 /3+2r(i + / 9) g ? )(fcfi)e b = b <*') = ~ 2*-pt( 2-p) H i-e {kR)ei{1 l3)t 

In the small-distance approximation Eqs. (12) and @) take the form (for —oo < x < L) 


v(x) = 


H(x) = ~ 


itan 7 ry(x) 

7r 

itan7rx(a;) 


In 


— OO 
2 


\x-t\ 


/z(t)dt 


d 2 

^2 J ln(|a; — t\)v(t)dt 


. kL 7 T 

In-h 7 H- 

2 1 2 i 


+ d y T(x) 


n(t)dt 


+ J c {x) 


(129) 

(130) 

(131) 


with x{x) defined in Eq. Q. 

In order that the contributions from large negative values of t will be small, the following asymptotics is required 


y,(t) ~ \t\ 71 , 

£—>■ — 00 


At) ~ ir 72 , 

£—>—oo 


7 j > 1 - 


(132) 
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Differentiating Eq. (951 and introducing functions (991 one gets the equations 


C+{x) = tan7rai 0(x)/+(x) — 


£_(x) = tan7rai 0(x)/_(x) + 


tan7rx(x) f L C +(t) 


£ 


£ 


—oo 
L 


t — X 


-d t 


tan 7 rx(a;) /' (~(t) 


t — x 


-df . 


For the incident field ( 128| ) 

/+(x) = \{d x + id y )'i>p(x, 0) = (1 - /3)bx~P, /_(x) = \{d x - id y )^p{x, 0) = pax 13 ' 

Introducing similar to Eq. 11021) analytic functions 


* ±(2) = si 


c±w 

t — z 


d t 


permits to find the general solution of Eqs. (133) and (134) (cf. Eqs. (106)) 

C\ sin7ro;i f L t a2 (L — t) ai 


$+(*) = 


z a2 (z — L) a i 2ni z“ 2 (z — L) ai 




'0 


C, 


sin 7rai 


$_C0 =_ 2 _i_ 

w z 1 ~°‘ 2 (z - 2tt\z 1 ~°‘ 2 (z - L) 1 - 011 


t — z 

r L t l - a *{L - t ) 1 ~“ 1 


/ 0 


t — z 


f-(t) dt 


(133) 

(134) 

(135) 

(136) 

(137) 

(138) 


Branches are fixed by requiring that fractional powers are real at real z > L. Imposing the correct behavior at large 
negative x (132), one concludes that for 0 < oq + a 2 < 1 (he. rj = 0) constant C\ = 0 and for 0 < a± + a 2 < 1 (i.e. 
77 = 1) constant Cq = 0. Notice that for opposite fluxes (i.e. when a 2 = 1 — a 1 ) the both constants are non-zero. 

The remaining integrals in Eqs. (137) and (138) reduce to the following ones 


t a (L — t) 


1 —a. 


Jo t- z 
For 77 = 0 one gets 


-df = 


sm na 


[z a (z — L) 1 a — 2 + (1 — a)L] , f — 

Jo 


t~ a (L — t) a , 7T 
v ; df = -— 


t — z 


sm 7 ra 


[z~ a (z - L) a - 1] . 


<F 


( 0 ) 


= (1 Q-i a 2 )b r ai ^ _ L)ai _ 

2i z a2 (z-L) a ' 1 y J 

1 - (ai+a 2 )o - r ai ( z _ - 0 + (1 - ai )L] . 

z 1 ~ a2 (z — L) 1_ai 2i z 1 ~ a2 (z - L) 1 -^ L J v 7 j 


(139) 

(140) 


^{Z) = -rj— 


For 77 = 1 


$W ( z ) = - Cl + | 2 ai —^ [z 1 " 01 (z - L) ai - z + ai L] , 

+ V > z a2 (z-L) a ^ 2iz a2 {z - L)“i L ’ J 


$ (1) (z) = - . ^ + , Q2 [z^-\z-Lf-^ - 1] . 

” ^ ' 2 lZ 1 - a 2 ( 2 ; - L)!-«i 1 ' ' J 


(141) 


Functions C±\t) are boundary jumps of these functions (cf. (103)). They have different forms depending on the cuts 

sin7rai F±\t), 0 <t < L 


c£ } (t) = 


sin7r(o!i + a 2 ) G±'(i), t<0 


(142) 


where 

• 77 = 0 


= ( L :° 1 ~ GfW = -(1 - a, - - f) "' - 1 


t “ 2 (L - f)“! 


( —t)“ 2 (i — t) ttl 


(143) 


Fi 0) (f) = - 


2iC 0 + (ag + a 2 )a[f - (1 - a\)L] q(o )_ 21(70 + (cti + a 2 )a[(-f) ai (^ ~ ^) 1 ai + t ~ (1 - cti)L] 

-/-I—rvo / T —rvi ’ ' ' ( J-M—rvo / T -l-\ 1—rv 1 


f 1 -“ 2 (L - t) 1 -"! 


(—t) 1_Q2 (L - t) 1 " 01 


16 







































• r) = 1 


n 1, w = -- 


- i)«l 


(i) ^ _ _ 2iCi — (2 — «! — a2)&[t — ai-L] ^(i) ^ _ 2iCi — (2 — — a 2 )^[( —0:1 (L — t) ai +t — aiL\ 

{i)n\ _ ( Q i +a 2 ~ l)a[{-t) ai -\L - ty~ ai - 1] 


= GL‘'(t) = - 


(—f)“ 2 (L - t) a ' 

\l-ai 


^"“ 2 ( 1 , - f) 1 -"! 
Functions /x(t) and i/(t) are 


(-t)!-“ 2 (L - 


(144) 


dd*) = ^(t) - CfW), ^ }, W = + C^W)- 


(145) 


For all values of 77 functions $±( 2 ) have singularities at 2 = 0 and 2 = 1 and decay as 2 7 with 7 > 1 at infinity. 
Therefore (which can also be checked by direct calculations) 


d^d t = 0, 


(146) 


and the only condition to fulfill is (as .F(O) = 0) 


/ J" (r) Mi)di = a 


The necessary integrals can be calculated by differentiation of the Euler integral and the answer is 

rL /UIX a _ 1 r(a 1 )r(a 2 ) 


ati+ct2~ 


J ( { °\t)ln(^pJdt = -ZL 
j C+\t)hi d t = —ZL 1 ~ ai ~°‘ 2 b 


r(«i + a 2 ) 

r(l-a 1 )r(l-a 2 ) 
r (2 Oil a 2 ) 


[2iCo + on (on + a 2 — l)a] > 


o:i+a2 —1, 


r(ai)r(a 2 ) 


£_Ah(M)d.- M ' -r (ai + 0> -ir 

L (t) d ‘ = r< r(2 -’a) - a 2 ) 2> [2iC ‘ - (1 - “' )(1 - 


where 


Z = 


7 r cos(7r(ai + a 2 )) sin 7rai 
sin(7r(a 1 + a 2 )) 


Condition (147) fixes values of constants C v 

2 iC 0 = ai(l - 07 - a 2 )oI + RbL 2(1 - ai ~ a2 \ 

2iC\ = -(1 - ai )(ai + a 2 - 1 )bL- TaL 2 (« 1 +“ 2 - 1 ), 

where 

F(1 - ai)r(l — a 2 )r(ai + a 2 ) 


i? = i?(ai, a 2 ) = 


T = i?(l - ai, 1 — a 2 ). 


(147) 


(148) 


r(ai)F(a 2 )r(l - ay - a 2 ) 

Calculating the derivatives one gets the following relations 
\{d x +\d v )^{x,y) = \{d x +id y )'bp(x,y) - 2i$+(f), \{d x - i3j,)^(a;,y) = |(d x - id y )^ 0 (x,y) + 2 i<F_( 2 ). 
From the above formulae it follows that in a small vicinity of the origin 

+ I4,)®‘»W) - ( 1 ~ - 2iC “ + ( “‘ + aMz - (1 “ Ql)i) 


2 a2 (2 - L) ai 


r l-a 2 ( 7 — 


(2 — L) 1- " 1 


(149) 


(150) 

(151) 


(152) 

(153) 

(154) 




(2-i) 1 
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From these expressions it is possible to calculate functions a,j and bj from the definition (151 and Eqs. (1501 and (|151) 


,(o) _ SiCo + ai(ai + a 2 )aL _ TU2 A T i- 2 u 1 -a 2 ,,(o) _ (1 - on - 02 )b 


= aL a2 H- L 

a i 


a = —e 1 ' 


aiL 1 -^ 

ll 2iC 0 - (1 - ai)(oi + o 2 )aL _ 


! — 
) U 1 ~ 


(1 - ai)L “ 2 


(155) 


02 -b 1 


= e 


oL“i - —L 1_2Q2_ai 
02 


1,(0)_iTTtt! (1 - an - a 2 )& 

’ 2 (l-a 2 )i“i ’ 


and 


(i) _ (ai + 02 ~ 1 )« ,(i) _ - 2 iCi + (1 - oi )(2 - oi - o 2 )& _ , T i- 


a = 


a\L 1 


> °i — 


„(1) _ piirau («1 + a 2 - l)a ^ _ 
“ 2 “ 6 ’ 2 “ 


„ = M/ 

(1 - ai)L Q2 

(i) _ _ i,rai 2iCi + ai(2 — ai — a 2 )b L 


Ta 


j^2ct\-\-cx.2 —2 


(1 - a 2 )L ai 


1 — ai 

—5L 1_ai 


(156) 


Ta 

1 02 


r Q;i+2q;2 —2 


Differentiating these expressions on Lj , taking into account that derivatives of a and b will contain additional smallness 
and, therefore, can be considered as constants, one finds 


• 7] = 0 


1 — a 2 


• T) = 1 


Pi d 

$1 Cl ) L —70 


P 2 ^2 

$2 C2 y c->o 


/?i £1 

di Ci ) L-t 0 


RL~ 


Q 2 (l - 02) x Q, — QC 2 — 1 



(157) 


(158) 


1-02 


02(1 — O 2 ) 


J^O'2—O'l 


— OL\ — 1 


OL 1 


L“ 2 ~ 


0 


02 


1 — Ol 


TLai+a2 -2 (1 - a2 )L“i-«2-i 


P2 e2 __ 

^2 C2 7 l-vo 


oi(l - ai) ^,.^! 


a 2 

Qi 

1-02 


TL a 1+a2 -2 _(]_ _ ai ) L <* 2 -on-l 


In all cases relation (37) is fulfilled with 


p = Ol^-Olj c 2 " °x 


(159) 


(160) 


(161) 


02(1 — 02 ) 

Coefficients <5i and 5 2 for i] = 0 and ei and e 2 for r) = 1 are zero in the leading order. Instead the explicit calculation 
of the next terms it is convenient to use Eqs. (411 by taking into account that 


(d x + id y )a = nb, n = ~k 2/3 2 1 2,3 e 17r/3 ^ + , (d x - id y )b = r 2 a, r 2 = 


(162) 


Comparing coefficients in front of a and b one gets the following limiting values 
• 77 = 0 

81 = — a2 (l ~ ~ r 2 e~ 17rai L~ ai ~ a2 , S 2 = — M 1 - Q - 1 - « 2 ) L - ai - a2 , 


and 


2(1 - oi)(oi + o 2 ) 

01(1 — 01—02) 2a 2 


9 2 — 7 


2(1 - oi)(oi + o 2 ) 


r 2 L 


9 3 = 


2(1 — o 2 )(oi + o 2 ) 

_ 2(1 ~ Ol) 2 °| 1 ~ Ol 


riL 2a2 , 


Ol 


1 — 0^1 — OL2 


(163) 


( 164 ) 


T _ 02 ( 1 - 01 - 02 ) T- 2 C-L 
n 2(1 — o 2 )(or + 02 ) 2 


/ 3 = — ^-RL~ 2ai 


02 


1 “ 2 -nL 2ai . 


1 — Qi — 0(2 
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V= 1 


_ (1 a 2 )(Ql + Ct 2 1 ) lwQl , ctl _)_ Q , 2 _2 

e l — 8 To \ rie ^ ! 

2ai(2 — ai — a 2 ) 


£2 = 


(1 ai)(ai + a 2 l) m i7rai rai+a 2 — 2 


-r!e ' 


and 


92 = - 

h = - 


2a-i 


1 — ai 
2 a 2 
1 — «2 


TL 2 “ 1-2 

rp^2a 2 -2 


ai 


ol i + a 2 — 1 
cc 2 

Q!i + 0:2 — 1 


r 2 L 2 ~ 2a \ 

r 2 L 2 ~ 2a \ 


2 o 2 (2 — aq — a 2 ) 

(1 - Oi)(ai + a 2 - 1) , 2a 2 


93 = 


a \(2 — ai — a 2 ) 


-nL 2 


(165) 


(166) 


/, = a-■»;)(“■+ °*-D ri£ 2 „ 1 - a . 

a 2 (2 — ai — a 2 ) 


In these calculations it has been taken into account that ( d x — id y )a and ( d x + i d y )b correspond to higher order terms 
in expansion (128), therefore they were put to zero in the leading order. 


The small-distance behavior of y = <5i e 2 follows directly from the above expressions. For 77 = 0 

M/ 


V 


L-> 0 


2 /££\ - 2 (a!+a 2 ) „ ; „2 


sin" tt(oi + 02) T (ct! + P 2 ) ci 7 r ( ai+ a 2 ) 

sin7TOi sin7ro 2 r 2 (ai)r 2 (a 2 ) 


(167) 


For rj = 1 the limiting behavior of y is given by the same formula but with substitution oq —> 1 — aq, a 2 —> 1 — a 2 . 


C. Two vortices at large distances 


The knowledge of the one-vortex solution (see Appendix [B]) permits to calculate the two-vortex case within the 
perturbation series when the distance between vortices is large. Consider functions A\(x) and Bi(x). In the lowest 
order, when the second vortex is absent, their asymptotics is determined by Eqs. (75) and ( |B15 ). The existence of 
the second vortex even at very large distance modifies these expressions to the following ones 


Ai(fi) —► \ —^—e ikr f[°\6i)e io ‘ iei+io ‘ 262 , B^x 1 ) —4 i/-^-e i * r ^ o) ( 0 1 )e iai9l+iaa ® 2 . (168) 

\xi | —voo V 7 ri kr |^i |—s-oo V mkr 

Here 9\ is the polar angle around the first vortex, X\ = (rcos^rsinUi), and 9 2 indicates the polar angle with the 
center in the second vortex, x 2 = (Rcos9 2 , Rsm9 2 ). It is assumed that vortices and cuts are such as indicated in 
Fig. 0 Functions ffj°\9 ) and 9 j°\ 8 ) with j = 1,2 are obtained from Eq. ( |B15| ) by extracting the factor e laS 


fj O) (0) = -ic 3 {aj)e lH+I7rQy/2 , s} UJ ( 0 ) = c 4 (oy)e" 

Notice that at the position of the first vortex 9 2 = 0 but at the second vortex 9\ = it. This choice of cuts has as a 
consequence that functions A 2 (x 2 ) and B 2 (x 2 ) with x 2 centered at the second vortex are given by a slightly different 
expressions 


„( 0 )/ 


^—inaj/2 


(169) 


A 2 (x 2 ) 


\x 2 1 —^00 V TTlkR 


(0 2 )e io ‘ l( ' 9l - n ' ,+io ‘ 2e2 


B 2 (x 2 ) 


\x 2 1 —>-oo V TTlkR 


2 ; e ifeii < 72 0) ( 02 )e iai _7r )+ io!202 . (170) 


The first order corrections correspond to re-scattering of these fields on the second vortex. When coordinates are 
calculated from the second vortex (x 2 = x\ — L) and L —>• 00, A\[x) has the following asymptotics 


lim AAx 2 ) 

L—t 00 


D(L)f[°\n)e' 


3 iai 0i +ia 2#2 g —ikx 


D{L) = 


AkL 


TTikL 


(171) 


According to Eq. (72) the scattering function for this incident field is T 2 ( 9 , 7 r) in Eq. (B17) without factor e‘“ 20 which 
is included in the above definition, 

•F 2 (0,7r) = F(0)e iaa ®, 


ismro ie/2 
1 ’ 2 cos( 0 / 2 ) 


(172) 


In (72) radius, R, is counted from the second vortex. To shift it to the first vortex requires to write R ss r + Lcos9. 
Therefore the full contribution to function Fi (0, L) in two lowest orders is (when r —► 00 9\ = 9 2 = 9) 

F\{0, L) = (/i ( ° } (0) + D(L)f[ 0) (ir)F 2 (9)e ikLcos0 ^) e i( “ 1+ “ 2)s . (173) 
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In a similar manner 


Gi (9, L) = (^ 0) (0) + D(L)g[ 0) (7r)F 2 (6>)e ifeLcos e ) e i(ai+Q2)e . 


(174) 


Derivatives of functions F\(Q, L) and G\{9, L) over L according to Eqs. (83) are (in the lowest order) linear combina¬ 
tions of F^\6) = f^(9)e^ ai+ol2 ^ e and 0^(9) = g^^e 1 ^ 1 '^ 02 ' 16 . Performing the calculations one finds that 


P 2 &2 

£2 C2 J L—¥ C 

For the second vortex 


/ £3(^1 ) c i v ( a ,- ao )/2 ic 3 (ai) c i„( a , +a ,)/ 2 ^ 

fcsin -noil 03(02) c 4 (a 2 ) 

2 _ 104(01) iw(Ql+a2 )/ 2 C / i(cx 1) a ^)/2 

\ C 3 (a 2 ) c 4 (a 2 ) j 


D(L)e v ‘ 


F 2 (9,L)= (/^ 0) (6»)+ D(L)/ 2 (0) (O)Gi(0)e _ifcicose ) e i ( ai+a3) ® _i,rai , 
G 2 (9,L) = (g^ {9) + D(L)g£ ] (0)G 1 (9)e~ ikLcose ^ _ 

Here Gi(0) indicates the scattering amplitude J-i(0, 0) with factor e lai6 removed 


Fi(8,0) = G 1 (d)e i “ ie , 


P m _ sm ™i —iS /2 
Gl{9) - 2sin(0/2) 6 ■ 


Differentiating them by L and using Eqs. (|83|) one finds 

fcsin 7rai 


Pi Si 

Cl Cl J o 


/ £3(0^2) j 7r ( Q , 1 —a ?)/2 ic 3 (£ 12 ) i 7 r( ai +q 7 )/2 \ 


C3(ai) 


c 4 («l) 


104(^2) -i^(q 1 + a ,)/ 2 £4(^2) i 7 r (q ,- a ,)/ 2 


\ c 3 (ai) 


c 4 (a 4 ) 


D(L)e~ 


These asymptotic values obey Eq. (371 with 


p = Qi(l - ai) c 2 ,ri ai 

a 2 (l — a 2 ) 


in agreement with (B14). 


The asymptotic behavior of variables (|50|) is 
V = S 2 ei 


L — 


— \k 2 sin7rai sin n a 2 .£/(£) 


z =/3 1 P 2 ;—> —^k 2 sm7raism7rot2D 2 (L) . 


(175) 


(176) 


(177) 


(178) 


L—> 00 


(179) 


(180) 


VII. CONCLUSION 


The main result of the paper is the demonstration that the problem of two Aharonov-Bohm vortices is integrable. 
As it is often in integrable systems, the exact solution is lengthy and tedious. 

The solution has been obtained by a generalization of the method used in Ref. j3] to solve scalar diffraction problem 
on scattering on a finite slit in 2 dimensions. The principal steps of the solution are the following. 

• Due to singular nature of AB interactions, wave functions are fixed uniquely by their behavior in small vicinity 
of the vortices and only a finite number of coefficients is necessary to reconstruct wave functions. 

• To find these coefficients it is useful to introduce singular functions Aj ( x ) and Bj ( x ) independent on incident 
fields with prescribed singularities at vortex j (see Eqs. (16) and ©)■ 


For the Helmholtz equation in the plane (and in other symmetric space as well) there exists a group of first 
order differential transformations which commute with the Lagrangian and cancels the incident field. 

When any of such transformations is applied to the exact wave function, the resulting function corresponds 
to zero incident field. But, in general, the transformed function becomes singular in one or many vortices. 
As all invariant operators are of the first order, these singularities can be compensated by a suitable linear 
combination of auxiliary singular functions Aj ( x ) and Bj (x ). In such a manner one gets a large set of equations 
which express certain derivatives of the Green function and the scattering amplitude through functions Aj ( x ) 


and Bj(x) (see Eqs. (21), (26), (33), (801). 
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• Specializing these relations to vicinity of vortex positions proves that certain derivatives of functions Aj (x ) and 
Bj(x) are linear combinations of the same functions (see Eqs. (24), (29), (30), (34)). 

• Coefficients in these relations are functions of vortex separations and by calculating commutators of different 
group transformations one obtains a system of non-linear equations for them (see Eqs. 


• All necessary coefficients can be calculated from a solution of the Painleve V equation (59) or (after a non-linear 
Backlund transformation) of the Painleve III equation (631. 


• As all equations are differential, to really use them it is necessary to know values of coefficients in a certain 
point. Analytically, one can calculate asymptotics of these coefficients in the limit L — > 0 (see Eqs. ( 120)-(122), 
(|l59l)-(fT60|),(fT63l>-([l67l)) and/or L ->oo (see Eqs. (fT75l) - dT80b ). 


The method of Ref. (4j used throughout the paper is quite general and flexible. Originally it has been used for 
solving certain integral equations, see Refs. [4] and [18]. As it is demonstrated in this paper, it can also be adapted 
to the problem of scattering on two AB vortices. Its generalizations for similar problems for the Klein-Gordon and 
Dirac operators in the Minkowski and Euclidean spaces (and probably in other symmetry spaces as well) seem to be 
possible. 


Notes added 


The principal ingredient of the above solution was the adaptation of the method of Ref. [3] to problems of singular 
AB vortices. After the paper has been practically finished, O. Lisovyy has remarked to the author that similar 
equations (even for an arbitrary number of vortices) had been derived by M. Sato, T. Miwa, and M. Jimbo in Ref. [9j 
in a different manner. That work is one in the long series of papers devoted to developments of the theory of holonomic 
quantum fields (see e.g. m and references therein). In Ref. [S], the authors constructed wave function with prescribed 
monodromy around a finite number of points. In two dimensions, monodromy transformations for the scalar equation 
reduce to the appearance of the phase factor e 2maj after encircling a point j which corresponds exactly to a AB 
flux line at this point. To get the necessary equations, the authors of Ref. [9] wrote the most general behavior of 
auxiliary functions Aj (x ) and Bj(x) in small vicinity of vortex positions as in Eqs. (38), (39). Computing the action 
of operators commuting with the Lagrangian and using the uniqueness as has been done in previous Sections one gets 
the same system of equations as above. 

The main differences of this paper and of Ref. S] is in the later the Euclidean space has been considered. So the 
Klein-Gordon equation (i.e. the Helmholtz equation (| 8 ]) with reversed sign of k 2 ) is 


(d 2 + d 2 ~k 2 )*(x) = 0 . 

The analogue of the radiation condition © in this case takes the form 

1 n -kR 
y/~R 


T(£) —► -A-e~ kH F{e). 


(181) 


(182) 


Therefore, complex conjugation of the solution does not change the correct asymptotics at infinity and wave function 
with all opposite fluxes (necessary in Section IV) is simply the complex conjugate of the initial wave function. In 
the Minkowski space used throughout the paper, complex conjugate turns out-going waves to in-going ones and is 
not an allowed transformation. For two vortices in Minkowski space these two functions are related by the inversion 
at the line connecting the vortices (see Eq. ( A12[ ) ). In general, wave functions with opposite fluxes appeared in the 
reciprocity relation should be calculated separately which effectively double the number of unknown variables. 

Another difference between this paper and Ref. ;9] is that in the latter the question of correct limiting values of 
necessary variables has not been discussed. Even for two vortices calculations of wave functions in the limit of small 
and large vortex separation is a complicated problem (see Section VI). For larger number of vortices it remains an 
open question. 

In general, interrelations of the theory of holonomic quantum fields and the AB problem is not widely known 
and fully understood (e.g. there is no reference in [9J to the paper of Alraronov and Bohrn, Ref. 0) and further 
investigation of this subject is of interest. 
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Appendix A: Uniqueness of solution and reciprocity relation for scattering on AB vortices 

The standard way of proving general statements about wave equation solutions is the use of current conservation. 
Let Ti and $2 be two solutions of the Helmholtz equation ([8]). The current conservation means that 

j) Jds = 0 (Al) 

where current J is 


J = ^ 2 d s ^i - ’Li^’La, 


(A2) 


and the integration is performed along any closed contour inside which there is no singularities of ’Ll, 2 - For two AB 
vortices possible contours of integration can be chosen as in Fig. |T] 

To cancel the current along the both sides of cuts due to phase jumps it is necessary to choose boundary jumps of 
'Ll and IL 2 differently. If 'Ll obey conditions ([ 2 ]) then to have zero current through the cuts function T 2 should obey 
the same conditions but with reversed signs of all fluxes 


’L 2 +(z,0)=e lma i'b 2 _( y x, 0), 3 y 'L 2 +(x ) 0) = e 2 maj d y ^ 2 -(x, 0) 


x £ Cj. 


(A3) 


When these conditions are fulfilled, it remains to check conservation of current along two other types of contours. The 
first consists of small circles around each vortex. If the both functions obey regularity condition ( fIo| ), the integral of 
the current over such circles tends to zero with decreasing of the radius. The last contour is the circle of large radius 
encircling all vortices (cf. Fig. [I]). Its treatment depends on the problem considered. 

If there exist two solutions corresponding to scattering on AB vortices then their difference, dT, obey the Helmholtz 
equation and all conditions p|)-(fTo|) with zero incident field. Choosing ’Ll = J'L and ^2 = <5\L* one concludes that 
the conservation of current implies that 


lim R 

R —> 00 


(<5’L*d fl <5'L - <5'L<9 i? <5’L)d0 = 0. 


(A4) 


But 


Idtfd’L - i/c<bL| 2 = Ids'Ll 2 + k 2 \S^\ 2 + ifc(c5’L*d fl <5'L - 5^d R 5^ . 
Using radiation condition ([ 9 ]) and the previous expression one concludes that if k 7 ^ 0 


(A5) 


lim R 

R —> 00 


|c5’L| 2 d^ = 0, 


lim R 

R —> 00 


|3ij5\L| 2 d^ = 0. 


(A6) 


Let v = )Th aj be the total flux of the vertices. Outside the circle of radius R which includes all vortices function J’L 
can be expanded in formal series on Hankel functions 


E A nH^ n (kr)e 


j.(v+n)<t> 


(A7) 


As Hankel functions decrease when r -> 00 as r -1 / 2 , from (A6| it follows that all |A„| 2 = 0. As these coefficients are 
coefficients of expansion over a complete set of functions one concludes that the only possibility is that <5\L = 0. In 
other words, the only solution obeying all conditions (|2|-([l0|) with zero incident field is identically zero. 


Similar arguments are used to find the reciprocity relation which relates the Green functions with interchanged 
positions of the source point and the observation point. For the scalar diffraction problems these two functions are 
equal but for scattering on AB-vortices one has to reverse all vortex fluxes. Let us denote the Green function for the 
scattering on vortices a = aq,..., a n by Gs(x, x') (0 < aj < 1). Then the reciprocity relation reads 

Gs(x,x') = G_s(x' ,x) (A8) 


where —a = 1 — aq,..., 1 — a n . 

The proof of this formula can be done as follows. By definition, the both functions obey the Helmholtz equation 


(A + k 2 )Gs{x,x') = 5(x — x'), (A + k 2 )G_$(x,x") = S(x — x"). 


(A9) 
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and on any c ut Cj they have opposite phase jumps. 
From Eqs. (A9| it follows that 


where 


G s {x",x')-G_ s {x',x") 



(A10) 


J(x) = G s {x,x')d s G_s(x,x") - G_ 5 (£,f ,, )9 ^ G 5 (x ) f , )) 


(All) 


with the integration is being taken over the same contour as above. As the both Green functions obey the same 
radiation conditions Q the current over a big circle tends to zero which proves reciprocity relation (A8). 

In [9] the Euclidean case, k 2 < 0, has been considered. As a consequence, = 'Ft(x), as in this case the 

complex conjugation does not change the asymptotic of wave function on infinity (’F ~ e~ kr ). For real k when T ~ e lkr 
complex conjugation contradicts the radiation condition <© and \F_ 3 (x) is not related immediately with tF^a;). For 
the problem of two vortices there exists an additional symmetry, namely the reflection in the line connecting two 
vortices. As such inversion interchange upper and lower parts of the cuts, conditions ([2| are now fulfilled but with 
opposite fluxes (in other words, fluxes are pseudo-scalars). Another method to check this relation is to consider 
small-a: behavior of wave functions Eqs. (10). It is clear that the inversion of y-coordinate is equivalent to reversing 
the flux, a —► 1 — a. Therefore up to a phase factor 

*a(£) = *-a(££) (A12) 

where the inversion S acts as follows 

S(x,y) = (x,-y). (A13) 


Combining it with Eq. (A8), one concludes that the following form of the reciprocity is valid for the two vortex 
problem 


G s (x,x') = G a (Sx',Sx). 


(A14) 


To determine constants t, in (36) one can proceed as follows. Close to a vortex with flux a situated at 0 the Green 
function behaves as in (151 


G(; x, x') —> a{x') {x + i y) a + b{x') (x — i y) 

x —>-0 


1 - 0 ! 


(A15) 


Assume for simplicity that the cut associated with this vortex lies on negative 2 :-axis. When the position of the vortex 
is shifted from 0 to SL along the cc-axis the new field close to the vortex remains practically unchanged 


Gsl(x,x') a(x') {x — SL + \y) a + b{x') (x — 6L — iy) 1 


>8L 


(A16) 


but a new portion of the cut from 0 to 6L appears. (Implicitly it is assumed that 6L > 0 so the length of the cut 
increases). From (A16) it follows that the difference of the field on the both side of such new cut is (for x £ [0, SL]) 


Gg^\x,x') — Gg L J (x,x') = 2isin7ra \a(x')v\(x) — b(x > )u 2 (x )\, 

J y=0 

d y Gg~^(x,x') — d y G^(x,x') =2is\m:a[a{x')iv' 1 (x)+ b{x')\u' 2 {x)] 

J y=0 


A-) 


where 


v\(x) = {SL — x) c 


V 2 {x) = {SL — x) 


1 -a 


(A17) 

(A18) 

(A19) 


These expressions determine the field on the both sides of a cut [0, SL], According to the Green theorem the field 
everywhere is given by the integral 


r&L 


G S l{x,x ) = G{x,x )+ / G{x,z) \dyG ( ^\x,x ) - d y G ( SL \x, x ) 
J 0 

r 8L 


J y=0 


dz 


d r G{x, z) G${x,x') - G k sl \x, x') 


(-)/ 


dz . 


J y-0 


(A20) 
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Here z denotes point (z, r) with r —> 0. 

Due to the reciprocity relation (A14) the behavior of the Green function G(x,x') when x' —>• 0 is the follows (cf. 

tm) 


G(x,x') —> a(Sx) (x 1 — \y') a + b(Sx) (x 1 + iy') 1 

x '^0 


It is plain that qj = 0 and 


27ra(l — a) rr 
ri=r 2 = - ; - oL. 


sin 7ra 


It means that one gets (361 with 

tj = 4a j (1 — a j). 

The same result follows from the one-vortex solution (see Appendix |b|) . 


(A21) 


;' + i y') 1 Q , x 1 = z, y' = r —> 0. 

Therefore in the leading order in 5L the Green function is 

G$l(x, x') ~ G(x , x') — 2 sin7ra a(Sx)[a(x')qi + b(x')ri] — 2 sin7ra b(Sx )[b(x')q 2 + a(x')r 2 \ (A22) 

where 

pSL pSL 

qi= [z a v[(z) + az a ~ 1 i / i(z)\dz, n = / [z a v' 2 (z) — az a ~ 1 V 2 (z)]dz, (A23) 

Jo Jo 

r*8L n8L 

<?2 = / \z 1 ~ a v' 2 {z) + (l - a)z~ a v 2 (z)\dz, r 2 = / [z 1_£ V((z) - (1 - a)z~ a v 1 (z)]dz. (A24) 

Jo Jo 


(A25) 


(A26) 


Appendix B: One-vortex solution and local constants 


The Green function for the AB problem with one vortex with flux a is 


G(£, x') 



n =—oo 


J\ n+a \(hr) e i ("+“) e H$ +al ( kR) e ~^+ a ^, r < R 
J\n+ a \(kR) e-^ + ^ H$ +a] ( kr ) e i (”+“) e , r > R 


Here x = ( 7 - cos d, r sin 9) and x' = (R cos (j),R sin <j>). 

When 0 < a < 1, terms with n = 0 and n = — 1 dominate in the limit \x | —> 0 and as |6j 




1 

r (l + J/ ) 


it follows that 


(Bl) 


(B2) 


Here 


with 


G(x,x') —> a(x')(x + iy) a + b(x')(x — iy) 1 “ 
| 5|->0 


i(x') = c 1 (a)Hi 1 \kR)e~ ia b(x') = c 2 (a)H[ 1 2 a (kR)e ~ i (“"W 


ci( a) = -- 


i k° 


02(a) = -- 


i k 


1 —a 


2 Q + 2 r(l + a) ’ ~ zv ~' 2 3 ~ a T(2 — a) 

From recursive relations for an arbitrary Bessel function Z v and the definition of the Hankel function [6] 


Z’ v (r) ± —Z v (r) = ±Z vT i(r), H™ = — -(J_„ - e~ im, J v ) , H™ = 

nr » l cin tt j/ x ' 


1 Sin 7 TV 


(B3) 


(B4) 


(B5) 


(B6) 
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it is straightforward to check that 

k sin 7 ra 

d L G(x,x ) =--- 


H^\kr)e' a6 + H^} a {kr)e~ Kl - a)e 


Eq. (21) states that the derivative over vortex position has the form 

dG(x, x') 


dL 


= a(x')A(x) + b(x')B(x). 


(B7) 


(B 8 ) 


As expected, for the one-vortex problem functions A(x) and B(x) are proportional to the Hankel functions 

A{x) = ca{a)H?} a (kr)S a - 1)e , B(x) = c 4 (a)H^ (kr)e ia9 , (B9) 

with 

03 (a) = — isin7raA: 1 ^ a 2 Q_1 r(l + a), 04 (a) = — isin7rafc“2 -Q T(2 — a). (BIO) 

The above expressions permit to calculate local constants for well-separated vortices. For the one-vortex solution one 
obtains 


(d x + i d v )A = g 2 B, 


<72 = ke' a7T —, (d x - \d y )B = g 3 A, 

C 4 


g 3 = — fc e _ia7r —. 

C 3 


(Bll) 


From these relation we conclude that diagonal terms g\ = 0 and <74 = 0. As these values are independent on L, they 
always remain zero. 

From Eqs. (45) one finds that 


mi = 1 — ai, 


Tii — —ct 1 , 


7712 = 1 - a 2 , 


n 2 — —a 2 . 


These constants are independent on L and it is these values that are used in the main text. 
According to the reciprocity relation (36) one should have 

A(x) = tb(Sx ), B(x)=ta(Sx). 


(B12) 


(B13) 


Here transformation S changes sigh of the second coordinate, S(x , y) = (x, —y). Its explicit form in polar coordinates 
depends on the choice of cut direction. S(r,9) = (r, 27 t£ — 9) where £ = 0 if the cut is along negative x-axis, i.e. 
— 7 r < 9 < 7r, and £ = 1 if the cut is chosen along positive x-axis, i.e. 0 < 9 < 2n. Comparing with the above formulas 
gives |f| = C 3 (a)/c 2 (a) = 47 ra(l — a) and 


t = 47 ra(l — a)e 


— 27ra:i£ 


(B14) 


From the asymptotics of Eqs. (B9) it follows that for the one-vortex solution functions F(9) and G(9) in Eqs. (75) 
and (80) are 


F(9) = c 3 (a)e i(Q - 1)e - i7r ( 1 -“)/ 2 , G(6) = c 4 (a)e i “ e - i7r “/ 2 . 


(B15) 


Using ( |80| ) and the above values of Cj(a) one finds that the scattering amplitude for the one-vortex problem has the 
form (it is assumed that the cut is along the negative x-axis, i.e. —n < 9,(f> < n) 


J-(M) = 


sm 7ra 


{(a—1/2)(6—(ft)-\-‘nra. sign cf> 


2 sin ((9 — (/))/ 2 ) 

In particular for cf> = ± 7 r the scattering amplitude is 

^($,±71") = — isin 7 ra 


gictO—iQ/2 

2 cos( 6 >/ 2 ) 


(B16) 


(B17) 


which agrees with ( 2 ] (with a correction nu). 
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